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Abstract: With the Pierre Auger Observatory we have the capability of detecting ultra-high energy neu- 
trinos by searching for very inclined showers with a significant electromagnetic component. In this work 
we discuss the discrimination power of the instrument for earth skimming tau neutrinos with ultra-high 
energies. Based on the data collected since January 2004 an upper limit to the diffuse flux of neutrinos at 
EeVenergies is presented and systematic uncertainties are discussed. 
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Introduction 

The nature and the production mechanism of the 
cosmic rays of ultra-high energy (UHE), above 
10^^ eV, is still unknown. All proposed mecha- 
nisms are expected to produce neutrinos. Classical 
acceleration process of charged particles in astro- 
physical objects create neutrinos through interac- 
tions with the radiation within the source region 
or with the Cosmic Microwave Background (GZK 
neutrinos) [1]. In other type scenarios they arise 
as direct or indirect products of supermassive par- 
ticles. The T neutrinos are heavily suppressed at 
production. In the scenario of neutrino flavor os- 
cillation and a maximal O23 mixing, the flavor bal- 
ance changes when neutrinos reach the earth. After 
travelling cosmological distances, approximately 
equal fluxes for each flavor are obtained [2]. Tau 
neutrinos that enter the earth just below the hori- 
zon, the so-called skimming neutrinos, may un- 
dergo a charged-current interaction to produce a r. 
When the interaction happens close to the surface a 
T can exit the earth and its decay in the atmosphere 
can produce an Extended Air Shower (EAS) de- 
tectable with the PieiTe Auger Observatory [3]. In 
the EeVrange, this channel has been shown to in- 
crease the prospect of detecting UHE neutrinos [4] . 



Search for neutrinos 

UHE particles interacting in the atmosphere give 
rise to EAS with the electromagnetic component 
reaching its maximal development after a depth of 
the order of 1000 g cm~^ and extinguishing grad- 
ually within the next 1000 g cm^^. After a couple 
of vertical atmospheric depths only the muons sur- 
vive. As a consequence very inclined showers in- 
duced by nuclei (or possibly photons) in the upper 
atmosphere reach the ground as a thin and flat front 
of hard muons. On the contrary, if a shower be- 
gins development deep in the atmosphere (a tau de- 
cay) its electromagnetic component can reach the 
ground and give a distinct broad signal. Therefore, 
the detection of very inclined showers with a sig- 
nificant electromagnetic component are a clear in- 
dication for UHE neutrinos. 

The signal in each station of the surface detector 
is digitised using FADCs, allowing us to unam- 
biguously distinguish the narrow signals from the 
broad ones and thus to discriminate stations with 
and without electromagnetic component (figure 1). 
We tag the stations for which the main segment of 
its FADC trace has 1 3 or more neighbour bins over 
the threshold of 0.2 VEM [5] and the area over 
peak ratio [6] is larger than 1 .4. The event is se- 
lected if the tagged stations fulfil the trigger condi- 
tion and they contain most of the signal. After this 
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Figure 1 : FADC traces from a station of two differ- 
ent real showers after subtraction of baseline and 
calibration. Top: moderately inclined (40°); bot- 
tom: quasi-horizontal (80°). 



selection, an almost pure sample of young showers 
is isolated. 

The next step uses the footprint of local stations 
included in the global trigger to select very in- 
clined showers. First a tensor is built using the sta- 
tion signals and the ground positions (in analogy 
to the inertia tensor) and the corresponding ma- 
jor and minor axes are used to define a "length" 
and a "width". Then, for each pair (i,j) of tanks, 
a "ground speed" is defined as di j/j At^.j |, where 
dij is the distance between them (projected onto 
the major axis) and |Ati.j| is the difference be- 
tween the start times of their signals. Horizontal 
showers have an elongated shaped (large value of 
length/width) and they have ground speeds tightly 
concentrated around the speed of light. In figure 2, 
we show the distributions of these discriminating 
variables for real events and simulated tau show- 
ers. The following cuts are applied: length/width 
> 5, average speed G (0.29,0.31) m ns^^ and 
r.m.s. (speed) < 0.08 m ns^^. We keep about 80 
% of the T showers that trigger the surface detec- 
tor The final sample is expected to be free of back- 
ground. 



Acceptance and neutrino limit 

Both the criteria to identify neutrino induced show- 
ers and the calculation of the ly^ acceptance are 
based on Monte Carlo techniques. The former uses 
the simulation of the shower development in the 
atmosphere as well as the detector response. The 
latter needs the simulation of the interactions that 
happen while the neutrino crosses the earth [7]. 

The total acceptance collected from January 2004 
until December 2006 with the Pierre Auger Obser- 
vatory is the time integration of the instantaneous 
aperture. 



Acc{E^) = / dEr 



AcCr{Er,hc) 



dt dxdy Ieff{Er,hc,x,y,Aconf 

JT J A 



where dN^ /dE^dhc is the flux of emerging rs and 
lef f the probability to identify a r. It depends on 
the energy of the t (Et), the altitude of the shower 
center defined 10 km after the decay point (he) 
[4], the instantaneous configuration of the detector 
(Aconf{t)), and the relative position of the shower 
footprint in the array (x, y). 

The Acc{Eu) is computed by Monte Carlo in two 
independent steps. First, the integral on time and 
area are performed using the simulations of the 
EAS and the detector, allowing us to account for 
the time evolution of the detector. The second 
step computes the integral on he and Er by adding 
Accr {Er , he) for all emerging r, given by the sim- 
ulation of the earth interactions. The statistical pre- 
cision due to the statistic of the Monte Carlo simu- 
lation is at a few percent level. 

The Monte Carlo simulations use several physi- 
cal magnitudes that have not been experimentally 
measured at the relevant energy range, namely the 
cross-section, the r energy losses and the r polar- 
isation. We estimate the uncertainty in the accep- 
tance due to the first two to be 15% and 40% re- 
spectively, based on Particle Distribution Function 
(PDF) uncertianties. The two polarizations give 
30% difference in acceptance. We take it as the 
corresponding uncertainty. The relevant range for 



30th International Cosmic Ray Conference 



0) 

> 
m 

o 10' 



SI 

E 



10 



V. 



10 



15 20 
length 



25 30 
/ width 




0.2 0.3 0.4 0.5 0.6 
Mean speed [m ns"^ ] 



i2 10' 



SI 

E 



10' 



10 



0.05 0.1 0.15 0.2 0.25 0.3 
RMS speed [m ns"^ ] 



Figure 2: Distribution of discriminating variables for neutrinos with an flux (histogram) and real events 
passing the "young shower" selection (points). Left: length/width ratio; middle: average of the speed 
between pairs of stations; right: rm.s. of the speeds. 



PDFs includes combinations of x and where no 
experimental data exist. Different extrapolations to 
low X and high would lead to a wide range of 
values for the v cross-section as well as the r en- 
ergy losses. The uncertainties on the low x regime 
as well as possible large ly cross-sections have not 
been added on the quoted systematics. 

We also took into account uncertainties coming 
from neglecting the actual topography around the 
site of the Pierre Auger Observatory (18%). We 
are confident on the simulations of the interactions 
undergoing in the earth at 5 % level. And we quote 
a 25 % systematic uncertainty due to Monte Carlo 
simulations of the EAS and the detector. 

Data from January 2004 until December 2006, 
which equate to about 1 year from the completed 
surface detector, have been analysed. In figure 
3, we show the collected acceptance on the anal- 
ysed period, for the most and least favourable sce- 
narios of the systematics. Over that period, there 
is not a single event that fulfils the selection cri- 
teria. Based on that, the Pierre Auger Observa- 
tory data can be used to put a limit for an in- 
jected spectrum K ■ with a known shape. 
For an E^^ incident spectrum of diffuse i^r, the 
90% CL hmit is E"^ ■ dN^r/dE^ < l.StEJj 
10~^ GeV cm~^ sr~^ s~^, where the uncertain- 
ties come from the systematics. The central value 
is computed using the v cross-section from [8], the 
energy losses from [9] and an uniform random dis- 
tribution for the tau polarisation. The bound is 
drawn for the energy range 2 10^^-5 10^^ eV over 
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Figure 3: The Pierre Auger Observatory accep- 
tance from January 2004 until December 2006 for 
the most and least favourable scenarios of the sys- 
tematics. 

which 90% of the events are expected. In figure 
4, we show the limit from the Pierre Auger Obser- 
vatory in the most pessimistic scenario of system- 
atic uncertainties. It improves by a factor ^3 in 
the most optimistic one. Limits from other experi- 
ments are also shown assuming a 1 : 1 : 1 balance 
among flavors due to the oscillations. 
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Figure 4: Limit at 90% C.L. to an E'^ dif- 
fuse flux of Vr at EeV energies from the Pierre 
Auger Observatory. Limits from other experiments 
[10, 11, 12, 13, 14, 15, 16] as wefl as fluxes for 
GZK V [17, 18] are also shown. For each experi- 
ment, the flavors to which is sensitive are stated. 

Summary and Prospects 

The dataset from January 2004 until December 
2006, collected by the Pierre Auger Observatory, 
is used to present upper limits on the diffuse in- 
cident Vr flux. The skimming technique is flavor 
sensitive and together with the configuration of the 
surface detector gives the best sensitivity around 
few EeV, which is the most relevant energy to ex- 
plore GZK neutrinos. The limit is still consider- 
ably higher than GZK neutrino predictions. Neu- 
trinos that interact in the atmosphere can also be 
distinguished from nucleon showers [7]. Hence, 
the Pierre Auger Observatory can explore UHE v& 
with two techniques that depend differently on v 
properties like flavour or cross-section. The Pierre 
Auger Observatory will keep taking data for about 
20 years over which the bound will improve by 
over an order of magnitude if no neutrino candi- 
date is found. 
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